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The field- and temperature-dependent de Haas-van Alphen oscillations spectrum is studied for 
an ideal two-dimensional compensated metal whose Fermi surface is made of a linear chain of 
successive orbits with electron and hole character, coupled by magnetic breakdown. We show that 
the first harmonics amplitude can be accurately evaluated on the basis of the Lifshits-Kosevich (LK) 
formula by considering a set of random walks on the orbit network, in agreement with the numerical 
resolution of semi- classical equations. Oppositely, the second harmonics amplitude does not follow 
the LK behavior and vanishes at a critical value of the field-to-temperature ratio which depends 
explicitly on the relative value between the hole and electron effective masses. 
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I. INTRODUCTION 

Due to their rather simple Fermi surface (FS), organic 
metals provide a powerful playground for the investi- 
gation of quantum oscillation physics. In that respect, 
the most well known example is provided by ^^-(BEDT- 
TTF)2Cu(NCS)2 which can be regarded as the experi- 
mental realization of the FS considered by Pippard in 
the early sixties for his model^. Such a FS, composed 
of closed hole orbits and quasi-one dimensional sheets 
coupled by magnetic breakdown (MB), yields quantum 
oscillations spectra with numerous frequency combina- 
tions that cannot be accounted for by the semiclassical 
model of Falicov and Stachowiak^. This phenomenon 
which is attributed, for this kind of FS, to either the for- 
mation of Landau bands or (and) the oscillation of the 
chemical potential in magnetic field has raised a great 
interestp^i^i^i^iii^i^. Nevertheless, the FS of numerous or- 
ganic metals is composed of compensated electron- and 
hole-type closed orbits^^, yielding many frequency com- 
binations as well, as far as Shubnikov-de Haas oscillations 
are concerne d^^i^^i^^ . In the case of a FS composed of two 
compensated orbits coupled to each other through MB 
but isolated from the other orbits outside the first Bril- 
louin zone (FBZ), it has been shown that the oscillations 
of the chemical potential can be strongly dampedi^ which 
could account for the absence of frequency combinations 
reported in the de Haas-van Alphen (dHvA) spectra of 
two-dimensional (2D) networks of compensated orbits in 
fields up to 28 T— . However, the FS considered in Ref. 
[13 which, to our knowledge, has no counterpart among 
the compounds synthesized up to now, does not provide 
a network of orbits and, therefore, do not yield Landau 
bands in magnetic field. 

The aim of this paper is to explore the field and tem- 
perature dependence of the dHvA oscillations spectra of 
an ideal 2D metal whose FS achieves a linear chain of 
successive electron- and hole-type compensated orbits. 
Such a topology can be relevant in compounds for which 




FIG. 1: (color online) Calculated Fermi surface of the Bech- 
gaard salt (TMTSF)2N03 in the temperature range below 
the anion ordering and above the spin density wave conden- 
sation, according to Ref. ITsl . Solid blue and red lines are 
compensated hole and electron orbits, respectively. 



the FS originates from an orbit with an area equal to 
that of the FBZ and coming close to the FBZ boundary 
along one direction, as it has been predicted for (BEDT- 
TTF)4NH4[Fe(C2 04)3}^. In this case, a large MB gap 
is observed at this point and a linear chain of successive 
electron-hole tubes separated with a smaller MB gap can 
be observed. Analogous topology is also realized in the 
Bechgaard salt (TMTSF)2N03 in the temperature range 
in-between the anion ordering temperature and the spin 
density wave condensationi^. We will focus on the field- 
and temperature-dependent amplitude of the first and 
second harmonics Ai and A2 of the oscillatory part of 
the magnetization Mqsc derived for such FS topology, 
which can be expanded as Mqsc = sin(27rzFo/6), 
where the fundamental frequency of the problem Fq and 
the dimensionless magnetic field b are specified in the 
next section. 
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II. MODEL 

We consider a 2D metal whose electronic structure 
consists of two parabolic bands with hole- and electron- 
character yielding a periodic array of compensated or- 
bits (see Fig. [2j). The bottom of the electron band 
is set at zero energy while the top of the hole band is 
at A > 0, with the possibility for the quasiparticle to 
tunnel through a gap between two successive orbits by 
MB. The total number of quasiparticles in the system 
is constant and, to lower the total energy at zero field 
and zero temperature, the area of the hole-band at the 
Fermi surface should be equal to the area of the electron- 
band which accounts for compensation. Indeed, if the 
hole band were completely filled, the zero-field energy 
would be higher than if one part of the quasiparticles 
were transferred from the hole-band to the electron-band. 
In the following, the effective masses linked to the two 
bands, m* and m^, can be different. As in Refs. 
the dimensionless field b and temperature t are given by 
b = B/B and t = T/T^ respectively, where B = h/eA^^ 
T= E/ks^ and Aq is the unit cell area. The effective 
masses and energies are expressed in free electron mass 
mo units and in units of E = 27rfi^moAo, respectively. 
Unit cell area of most organic metals is in the range 100- 
200 yielding B and T values of few thousands of Tesla 
and Kelvin, respectively. Therefore, realistic experimen- 
tal conditions yield small values of b and t compared to B 
and T, respectively. On the contrary, the ratio which 
is the relevant external parameter for perfect crystals, is 
given by b/t = {eh/mokB) B/T. Its value is close to 
the B/T ratio achieved in experiments since eh/mokB — 
1.34 KT~^. Given an energy E^ the areas of the electron- 
and hole- type surfaces are respectively = 27rm*£^ 
and Sh = 27rm^(A — E)^ which are both quantized for 
closed orbits in the Brillouin zone. The zero field Fermi 
energy /io is given by the condition of compensation 
5'e(/io) = Shil^o) hence /io = [^h/i^t + ^Dl^- The 
fundamental frequency of this system is therefore equal 
to Fo = Se{fio)/27r = Sh{fio)/27r = m>* A/(m: + m*). 
Each Landau level (LL) has a degeneracy b per sample 
area. The spectrum of such chain of coupled electron-hole 
orbits and the quantization of the energy are determined 
by semi-classical and conservation rules of the wave- 
function amplitudes at the junctions where the quasipar- 
ticles tunnel and across the boundaries of the Brillouin 
zone. In particular, the amplitude of the wave- function 
at different points of the Fermi surface (see Fig. [2] for 
notations) satisfies the following rules: its phase is pro- 
portional to the area swept by the quasi-particle around 
the trajectory divided by b. We will note cTg = Se/2b 
and ah = Sh/'^b the phases around half the electron and 
hole orbits, respectively. Also, at every junction there 
is a probability ip for the quasi-particle to tunnel and a 
probability q to be refiected. Finally we add a Maslov in- 
dex i = at the vertical extrema of the orbits (cross 
symbols on Fig. [2]). With these rules, we can write the 
relations between the wave-function amplitudes a's and 





FIG. 2: One dimensional chain of coupled electron and hole 
orbits. The Brillouin zone is delimited by the dotted lines. 
Amplitudes are chosen at specific points on the trajectory 
(see text). 



ai = iex.-p{i(Je){ipcy -\- qP)^ 

Oi2 = ipai + qP2, OL = i exp(-icr/i)a2 

(5x = qai + ipp2, (1) 

(32 = iex.p{-iah){qa' ^ ip(3'). 

As boundary conditions, we impose that the ampli- 
tudes across one Brillouin zone are identical up to an 
arbitrary phase 6: ol = aex.p{i6) and = f3ex.-p{i0). 
Solving this system of closed equations, we obtain the 
quantization of the energy which satisfies the spectrum 
equation 



1 -\- q^ exp(2i(7e) — p'^ exp i{ae — ah -\- 0) 
1 -\- q^ exp(— 2i(j/i) — exp i{ae — ah — ( 
-p^q^ exp i{ae — ah) + exp(— 2zcr/i + iO) 
expi((Je — ah) + exp(2zcre — i9) (2) 



We can distinguish two limiting behaviors. For q = 
1 (or, equivalently, p = 0), the spectrum reduces to 
[l+exp(2z<Je)] X [l+exp(— 2z<j/i)] = which corresponds to 
the discrete LL of two independent electron and hole or- 
bits: Se{En) = 27r6(n + l/2) d^nd Sh (En) = 27r6(n + l/2), 
with n positive integer. In the opposite case where q 
= 0, the spectrum corresponds to non-quantized phases 
ae{E) — ah{E) = ±6 since no closed orbit exists in such 
case (open chain) and therefore magnetic oscillations van- 
ish, leaving the place to a continuum of states. To study 
numerically the spectrum given by Eq. ([2]) , we first deter- 
mine the periodicity in energy of the discrete LLs. For 
that purpose, we assume that mJ^/m*=(7o/po where qo 
and po are coprime integers that can be as large as needed 
to approximate the ratio of the two masses. It is indeed 
useful to express this quantity with a rational number 
since it allows to consider as necessary only a finite set of 
solutions of Eq. ([2|). The minimal periodicity Te of the 
spectrum is, in this case and for a given phase 6>, equal 
to Te = 2bpo/ml, which is twice the periodicity calcu- 
lated for an isolated system made of one hole and one 
electron band^^. The number of solutions inside each in- 
terval of width Te is found to be equal to 2(po + ^o) (this 
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number is conserved when p varies and can be counted 
exactly for p = 0). Given a set of solutions Eeh{q^n^O)^ 
n = 0, • • • , 2(po + ^o) — 1, we introduce a cutoff function 
(Pc{E) such as (pc{E) = 1 for E larger than a characteris- 
tic energy Ec and equal to exp[—c{E — Ec)'^^] for E < Ec, 
where 6 is any positive integer greater than 1 (we take 
(5 = 4 in the simulations which gives a very smooth cut- 
off function), and c a positive parameter determined self- 
consistently. This function has the property of preserving 
the physical features near the Fermi surface and the GS 
energy is in particular finite since now the corresponding 
hole spectrum is bounded for large and negative energies. 
The LL density of states pc{E) takes the following form 



p27T 2(po+Qo)-l oo 

n—O k— — c<. 



b r"" 

^c{E)5{E - Eeh{q. n, 9) - RTe) (3) 



Given Ec^ the positive parameter c is found to be solu- 
tion of the equation of conservation at zero temperature: 



Neh = 



dEpe{E) 



^""^-oo Jo k=-oc 

ME)S{E - Eeh{q, n, 0) - kTs) (4) 



where N^h is, as mentioned before, the total number 
of quasiparticles in the Canonical Ensemble. In the nu- 
merical simulations, /io(^) is taken as the first LL lo- 
cated below /io- We choose A^e/i, which is arbitrary, as 
a multiple of the characteristic zero field energy den- 
sity (m* + Tnl^)fio (which is proportional to the areas of 
the Fermi surface). We will also choose in the following 
Ec = -2 and Neh = S{ml + m^)/io, and 10 values of 
for the integral evaluations. For each value of the field 
6, the parameter c defined from Eq. (|4]) is unique and 
determined self-consistently. Then we can compute for 
example the GS energy AEq 



E^c{E)S{E - Eeh{q. n, 0) - kTs) (5) 
and the free energy AF 



AF 



"2^ 



/oo />27r po+go-l 00 

dE dO ^^(^) 
J^ n=0 

log[l + exp;5(/i - E)]5{E - Eeh{q, n, 0) - kTE) 

+ NehP- (6) 




FIG. 3: (color online) Fourier spectra of the oscillatory mag- 
netization (displayed in the inset) for various q values. The 
effective masses are me = 1 and rrih = 2.5. Fq is the funda- 
mental frequency (see text). 
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FIG. 4: (color online) b/t dependence of the chemical poten- 
tial normalized to its value /io in zero- field for various q values. 
The effective masses are me = 1 and ruh — 2.5. 



The chemical potential ji = ii{t^ h) is calculated from 
Eq. by extremizing the free energy dAF/dja = 0. The 
magnetization Mqsc = —dAF/db = x^dAF/dx, where 
X = 1/6, has to be independent of the parameter Ec for 
Ec far away from the chemical potential or at energies 
large compared to the Landau gap for a given magnetic 
field. We have checked, for different values of Ec with 
A = 1, for example Ec = —1,-2,— 4, and for a large 
range of fields, that the resulting magnetization is stable 
in this procedure. 

Examples of field-dependent dHvA and chemical po- 
tential oscillations deduced from the numerical resolution 
of Eq. [6] are given in Figs. [3] and [H respectively. These 
data are derived assuming a constant tunnelling probabil- 
ity p = y^L^^q^ even though p is field-dependent. Nev- 
ertheless, a qualitative conclusion can be derived from 
these data which otherwise will be used in the following 
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for the study of the Fourier amplitudes: the oscillations 
of the chemical potential are weak, owing to the large h/t 
range considered, and decrease as q decreases. This point 
which is in agreement with the conclusion derived for iso- 
lated orbits^^ confirms that the field-dependent chemical 
potential oscillations are strongly damped in 2D compen- 
sated metals. 
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III. AMPLITUDE OF THE FIRST HARMONICS 
IN THE LK APPROXIMATION USING 
RANDOM WALK ANALOGY 

In this section, we compute the expression for the 
first harmonics amplitude ^, corresponding to the fre- 
quency Fo, within the LK approximation. We will check 
numerically that this approximation fits very well the am- 
plitude Ai extracted from the spectrum computed in the 
previous section. We first notice that there are an infi- 
nite number of orbits along the chain contributing to this 
amplitude (see Ref. [13 ) since a semi-classical trajectory 
around a successive electron and hole pockets gives a zero 
frequency contribution. Therefore these zero frequency 
orbits can in principle be added to any orbit which has 
a Fo contribution to form new Fo-orbits by juxtaposi- 
tion of continuous trajectories. These orbits yielding the 
frequency Fq can be classified by their successive masses 
me(l) Iml + (/ - 1)^^ or mh(l) = (I - l)^e + 
where / is a positive integer. The contribution of orbits 
with large effective mass is negligible at low field or high 
temperature {h/t small) but have a significant value in the 
large h/t limit where all the orbits have to be taken into 
account. Indeed the field- and temperature-dependent 
damping factors which are defined for a given set of ef- 
fective masses ^*(/^) by the formula 



sinh(27r2m*^^^t/6) 
are close to unity in this limit. Dingle factors 



(7) 



RD{ml^^ytl^^^) = exp(-27r^m*(^)t*(^)/6) (8) 

where = moAo/4:7r'^hre[h) are the reduced Din- 

gle temperatures, can also be added to account for real 
crystals in which the relaxation times r^i^h) have finite 
values. In the following we will assume first that these 
two relaxation times are negligibly small to simplify the 
calculations. The existence of an infinite set of orbits 
contributing to the first harmonics leads us to define a 
closed random walk on the chain in the Brillouin zone 
{xi}i=o,2n, with origin and end at xq and X2n, respec- 
tively, with the condition xq = X2n = 0. These coordi- 
nates take integer values (either negative or positive) and 
define precisely the pocket inside which the quasiparticle 
is located. We chose Xi with i even to be the positions 
of the electron bands and i odd for the hole bands. A 



FIG. 5: Example of a trajectory contributing to the first am- 
plitude, with fundamental frequency Fo and mass me(l) = 
2m* +771^. The successive coordinates {xi}i=o,..4 of the quasi- 
particle are given below the figure. Here = (see text). 
a -I = —1 since the particule is going backward when it is 
reaching the location xq = X4 from x-i = X3 on a periodic 
orbit. 



closed path has an even number of steps 2n. For a given 
Xi, the particle can also orbit a number > of times 
around the surface before going to the next band. An 
example of a typical closed trajectory is given in Fig. [5l 

We can rewrite coordinates Xi by the mean of for- 
ward/backward variables cr^ = ±1 such as Xj^ — z — X I 
ai{xi-i — Xi-2)' Here cr^ = 1 when the particle is going 
forward on its path (in the same direction) and = — 1 
when it is going backward in the reverse direction (see 
Fig. [5] for an explicit set of a^'s). The path is moreover 
made periodic by imposing Xi = Xi±2n (and the same for 
variables (Ji). On Fig. [5l X3 = x_i and the particle is 
moving on the reverse direction when it is reaching the 
point Xq = X4 along the portion of path (x_i,xo,xi), 
which implies that ai = —1. 

It will be also useful in the following to introduce new 
periodic variables yi = Xi — xi-i = ±1 which satisfy the 
simple relations = (Jiyi-i and cr^ = yiyi-i. Let us 
now calculate the total number of possible orbits with 
a given effective mass me{l) (or equivalently mh{l) for a 
hole orbit since the problem is symmetric by exchange 
of electron and hole masses) and their contribution to 
the amplitude of the first harmonic A^^ with the fun- 
damental frequency Fq. If the particle is going forward, 
there is an amplitude equal to ip to be transmitted from 
one band to the next one or ipq if the particle is going 
backward (it has first to perform exactly one reflection 
with amplitude q on the band edge before being trans- 
mitted back through the previous junction). In term of 
variables (7^, this is equivalent to write this partial am- 
plitude as ip[{l + (Ji) + (1 — (Ji)q]/2 = ip^e^^\ with 
= Then the LK amplitude corresponding to 

all electron-like contributions can be written as 



in: ^ me{l) 



2(i-l) 



^ 1 ^ (_p2)«q«+2E?:o'%e^^-i'"' 



n=l {ni,ai} 
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The sum over / corresponds to all possible masses, and 
the sum over n > 1 the number 2n of possible tunnellings. 
The first term / = 1 in Eq. (|9]) is the simplest orbit of 
the expansion: a closed trajectory around one electron 
pocket with MB damping factor and effective mass 
m*. The sum over n is limited to 2(/ — 1) which corre- 
sponds to the extremal trajectory, i.e the particle of mass 
me{l) visits a maximum of (/ — 1) successive electron and 
hole pockets outside the first initial electron pocket before 
going back, / — 1 being also the maximal linear extension 
of the path. This implies that all the rii are zero for this 
case (see Fig. [5|). The factor 1/n takes into account of 
the orbit symmetry by circular permutation of the co- 
ordinates {xi}. Finally, the sum over the set {rii^ai} is 
constrained by the boundary conditions xq = X2n and by 
the facts that the frequency is set to Fq and the effec- 
tive mass is /m* + (/ — This implies the following 
conditions on {rii^cri}: 



n-l 



n—1 

n 1 - o-2j+i 



j=0 j=0 

L 

j=o j=l 



n-l 



i-i = 



E 



1 - 0"2j 



(10) 



(11) 



The constrained sum can be transformed using two 
Kronecker integrals around the complex unit circle 



l>2 ' 



(12) 



2(/-l) 



{ni,ai} 



The sums over the {n^ = 0, . . . , oo} can be performed 
exactly, since q is less than unity. We obtain 



ml 



l>2 



me(l) 



2(Z-1) 



n=l 



2i7rz J 2i7rz' 



-l-\-3n/4 



I ^/-/+l+3n/4 -|_ 



(1 - ^E;roV2, + i/4 (1 - g2^/)n ^/E7=i^2,/4 

The condition for a closed path is given by Y^'j^i Vj = 
0, with the periodic conditions yo = which imposes 
in the previous sum the introduction of another Kro- 
necker function. The previous quantities depending on 
= yiUi-i can be rewritten in term of the ^^'s only. 



Using the fact that = the sum over the terms 

involving these variables is given by the quantity 



^^e-E?=i2/. X 



27r 



^0= E 

{y^=±l}' 

n-l 

Yl {zq^)-y^oy2j+i/^^^'q2yy2j+iy2j+2/^ (^4) 

This is the partition function for a one dimensional 
periodic Ising model with an imaginary field and al- 
ternate complex coupling. Setting the 2x2 trans- 
fer matrices Pyy^ = ^zq'^)-yy' /^e'^'^y+y'^^'^ and P^^, = 
^^fq2yyy'/4^iaiyW)/2^ Zq cau be Written as a trace over 
the product of operators {PP')^: 



(15) 



where \± are the eigenvalues of the matrix PP' . The 
amplitude can then be written as 



^.K(^) ^ g^i?(mg) ^ fi(m (0) ^ 
ml ^ me(0 



2i7iz T 2i7iz' 



-l+d,n/A /-^+l+3n/4 



(l-q^zY (l-q^z'Y^^' 

It is useful to express complex vectors z = e*^ and z' = 
e*^ by their angles 9 and 0' . Also, setting = y^e*^/^, 

= y^e*^ C = 1/C and C,' = l/C^ we can express the 
eigenvalues as 



A± 



CC'cos(2a) + CC', 

-(C'-C')(C''-C'') 

-(l-g2^)(l-gV)/(g2v^). 



(17) 



Then 



n/2 i 

^o = 2EEC'nC>"-V (18) 

i=0 j=0 
n/2 i n-2j 

= 2EE E ^^rC^Ct2,(CC')"-'^-'S^cos'=(2a) 

2 = j = k = 

where we introduced binomial coefficients = 
n\/k\{n — k)\. In the last line, we can perform the in- 
tegral over a of the integrand cos^(2a). The resulting 
integral is not zero only for k even: 
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L 



da 
2^ 



(19) 



TABLE I: First values of coefficients S{l,n) representing the 
number of non- equivalent orbits for a given mass me(h) (0 with 
2n magnetic breakdowns, 1 < n < 2(/ — 1). 



We finally obtain 



n/2 i n/2-j 
i=0 j=0 /c=0 



22k 



The amplitude A^^ can be rearranged using the pre- 
vious results like 



2{l-l) 



l>2 

n/2 i n/2-j 



rrieil) 



(20) 



E 2^^^^ EE E c^cic^^p^ 



i=0 j=0 /c=0 



dz f dz' z' 

X 



2k 

l-\-n—j — k y/— Z+l+n— j — fc 



X 



2fe 22/c 

2n— 4j— 4/c 



2i7rz J 2z7rz' (1 - q^z^-^ (1 - q^z'Y'^ 



To compute the last two complex integrals, we use the 
relation for any positive integers (a, /?) 



dz 



2mz (1 - q^zY 



2a 



(21) 



then 
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(22) 



where we defined the following combinatorial quanti- 
ties 
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FIG. 6: (color online) b/t dependence of the first harmonics 
amplitude. Solid symbols are deduced from the numerical res- 
olution of Eq. (|6)). Solid lines correspond to the LK approx- 
imation given by Eqs. (|22|) to (|24|) . Dashed lines correspond 
to the first order term (/=1 in Eqs. (|22|) ). The effective 
masses are me=l and m/i=2.5. 



i/2 



a/2- 



^(^'-) = ^EEE 



i=0 j=0 /c=0 



22k 



r^2ir^jr^2k r^k ^l — n-\-j-\-k^l — n-\-j-\-k—l 
^n^i^n-2j^2k^l+k-l ^l+k-2 



(23) 



These positive integers S{l^n) count the number of 
non-equivalent orbits for a given mass me(^)(/) and for 
which the quasiparticle is visiting 2n successive pockets. 
In table I we have reported, for information, the first 
numbers for increasing values of / from 2 up to 7, us- 
ing relation ([23|) . For a given /, n is taken from 1 to its 
maximum value 2(1 — 1). 

Finally, the total amplitude for the first harmonics is 
given by the sum of hole and electron contributions: 



TT 



Af'^(e)+^t^(/i) 



LKi 



(24) 



The first harmonics amplitude is determined by the 
set of Eqs. ([22]) to ([24|) . Examples are given in Fig. 
[6] (solid lines) and compared to the numerical resolu- 
tion of the Fourier spectrum of the magnetization taken 
from Eq. (|6]) (solid symbols) and for different q val- 
ues. Since the q is fixed, these data can be regarded 
as temperature-dependent amplitudes at a given field. 
An excellent agreement between numerical results and 
Eqs. ([22]) to ([2^ is observed, even at high magnetic 
field. Dashed lines in this figure are the contributions 
of the first order terms (Eq. [231 reduces to A^^ = 
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q'^Fo[R{ml)/ml + i?(m^)/m^]/7r within this approxima- 
tion). These terms, which only takes into account the 
basic orbits with the lowest effective masses (mj and mf^) 
are strongly dominant since only a small difference (less 
than 10 %) is observed in the high b/t range. However, as 
discussed below, the high-order terms with higher effec- 
tive masses can have a significant influence on the evalu- 
ation of the effective masses from experimental data. 

In the case of real experimental data collected on q- 
2D compensated metals, an effective mass m* is deduced 
from the temperature dependence of the amplitude of 
the first harmonics Ai at a fixed magnetic field. In such 
a case, it is implicitly assumed that either the effective 
mass of electron and hole orbits is the same or, oppo- 
sitely, that only one orbit contributes to the considered 
Fourier component because the other has a much larger 
effective mass. In addition, the contribution of high or- 
der orbits (/ >1) is neglected. Taking into account the 
Dingle damping factor (see Eq. [8]), we can rewrite the 
LK formula as: 
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FIG. 7: (color online) b/t dependence of the effective mass 
deduced from Eq. (|26p with the parameters relevant to data 
in Fig. El 



y (X Rmb — 



*t*/b 



sinh(27r2m*t/6) 



(25) 



where y = Aib/t, Rmb is the relevant MB damp- 
ing factor and t* is the reduced Dingle temperature. 
If the magnetic field is fixed, Eq. ([25]) reduces to 
y (X 1/ sinh(27r^m*t/6). Therefore, the effective mass can 
be extracted by considering the following combination of 
derivatives: 



(27r2m*)2 



1 % l2 1 



ydit/b)i ydHt/b) 



(26) 



Figure [3 displays the b/t dependence of m* correspond- 
ing to the data in Fig. [6] which stands for a perfect crys- 
tal (t*=0). For a given q value, the results are scaled as 
m*/min(m*, m^). For q=l^ data analysis based on Eq. 



yields m*/min(m*, m^)=l and Y^mjm^/min(m* , 
m^) in the low and high b/t limit, respectively. This 
point is further supported by considering the mass plot 
in Fig. [8] which yields a straight line at high t/b. 

For g < 1, a strongly non-monotonic b/t dependence 
of m* is observed in Fig. [71 This behavior is linked 
to the high order terms (/ > 1 in Eqs. ([22|), ([23]) ). 
The g-dependence of this behavior is also strongly non- 
monotonic since numerous zeroes can arise in the coef- 
ficients involved in Eq. ([22]) as q varies. Nevertheless, 
the effective mass variations are damped for real crystals 
with finite Dingle temperature, as demonstrated in Fig. 
|9l for (7=0.6 where it is assumed that = for sim- 
plicity. Assuming a fixed magnetic field of 30 T (which 
yields a MB field Bo=30.6 T for (7=0.6), r/6=0.01 stand 
for a good crystal with a Dingle temperature Td=0.4 
K. Oppositely, t*/6=l, for which m* is always close to 
min(m*, m^) corresponds to an extremely bad crystal for 
which uJcT = 6/27rm*t* (where the cyclotron frequency 
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FIG. 8: (color online) t/b dependence of y = Aib/t deduced 
from Eqs. {|22]) to (|24)) for m*=l, m^=2.5 and q=l (solid 
symbols) and best fits of Eq. (|26)) to this data in the low 
(blue solid line) and high (green solid line) t/b range. The 
inset displays the low t/b range. 



cOc = 27rhb/mom^Ao in our units) is much lower than 1 
at experimentally accessible fields. 



IV. PROPERTIES OF THE SECOND 
HARMONICS IN THE CANONICAL ENSEMBLE 

In this section, we study the t/b dependence of the 
second harmonics. If the analytical expression is difficult 
to obtain for general parameter (7, it is however possible 
to obtain some information in the limiting case (7=1, 
which has already been considered in a previous studyi^. 
In particular, we have seen that, in this limit, the free 
energy AF, for a compensated metal, is given by the 
difference between the Grand potentials of the electron 
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FIG. 9: (color online) b/t dependence of the effective mass 
deduced from Eq. ()26p for me = l, m^=2.5, ^=0.6 and various 
values of the Dingle temperature (t*) indicated on the curves. 



and hole bands, respectively: 



2 CO 
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(27) 
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exp(2i7rnFox + 2i7rnWeG{x)) 
exp(2z7rnFox — 2i7rnWhG{x) 



with 1(7, 



'e(/i) 



e(/i) 



/(m* + m^). We make the fur- 
ther approximation in the exponential parts of (|3Q|) , 
and which is valid at low temperature, that G{x) can 
be truncated to the first term G{x) ^ [i?(m*) — 
R{ml)] sin(27rFox)/7r so that 



CO 

E 



J,n(^ci^e(^))exp(2i7rmFo^) (31) 



where ae(/^) = 2me(/^) [i^(m*) -i?(m*)]/(m* +m*), and 
are the Bessel functions of integer order. Putting ex- 
pression ([3T]) inside ([3Q|) . we then select, in order to iso- 
late the second harmonics, integers such that n + m = ±2 
and n — m = ±2 for the electron and hole contri- 
butions, respectively. Expanding the magnetization in 
terms of Fourier components Mqsc = sin(27rFo) + 
A2 sin(47rFo) + . . . , we find that the coefficient A2 satisfies 
the relation 



The chemical potential /i satisfies the self-consistent 
equation dAF/dfi = given by 



Mo 



- ml ^ irn 



R{nml) sin(27rnm 



We follow the appendix B of reference (^l3) to extract 
(for small field values) the second harmonics from equa- 
tions ([27t [28|) . The oscillatory part of the magnetization 
is given by Mqsc = —dAF/db = x'^dAF/dx (we remind 
that X = 1/b). We then introduce the periodic function 



^(^) = Yl —l^i^^h) - R{nml)] sin(27rnFox)(29) 

n=l 



SO that the chemical potential Eq. ([28|) can be ex- 
pressed as /i = /io + bG{x)/{ml + ^/^), and the magneti- 
zation 



7TA2 _ (R(ml) - R{ml)f 



2Fo 2(m* + ml) 

[R{{n + 2)ml) - R{{n + 2)m*)][i?(nm^) - R{nm, 



(32) 



7 Y.^^[R{nml)sm{2^nml^) 
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n>l 



n(n + 2){m% + m^) 



El \R{nml) -r / X -r / XX 
— ZT^ — \Jn-2{nae) - Jn+2{nae)) 



L ml 



(28) 



1 — [Jn-2[nah) - Jn+2[nah)) 



In Fig. [To] is plotted the expression ([32]) together with 
the numerical results obtained by solving directly the 
spectrum ([2]) and extracting the second harmonics am- 
plitude for different values of q (see section [2]). The two 
results at g = 1 agree quite well for a large domain oit/b. 
Remarkably, the numerical data show that, whatever the 
q value is, the amplitude vanishes at a unique value oit/b 
depending on the ratio of the two effective masses only. 
In the example given by m* = 1 and m^ = 5/2, the am- 
plitude vanishes for t/b 0.1. This behavior is strikingly 
different from that of Ai = ^, considered in previous 
section, for which no singularity is observed. This reflects 
the fact that the deviation from the LK approximation 
appears only in the second harmonic. 
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FIG. 10: (color online) t/b dependence of the second har- 
monics amplitude A2 for various q values. Solid symbols and 
solid line correspond to the numerical resolution of Eq. [6] and 
to the solution of approximation (|32)) for ^ = 1, respectively. 
Here m* = 1 and = 5/2. In the insert is plotted the crit- 
ical value of t/b at which A2 vanishes as function of the ratio 
ml/ml. 



V. SUMMARY AND CONCLUSION 



The spectrum for one-dimensional chain of compen- 
sated orbits has been calculated. As it is the case for 
two isolated orbits, the field-dependent oscillations of the 
chemical potential are weak. In turns, the LK formal- 
ism can be applied, provided MB orbits, which in such 
a network contributes to the fundamental Fourier com- 
ponent amplitude, although with higher effective masses, 
are taken into account. The resulting high order terms 
(/ > 1 in Eqs. ([22|) and ([23]) ) can lead to apparent 
temperature-dependent effective mass for clean crystals 
in the high b/t limit in the case where only one effective 
mass is considered for the data analysis, as it is usually 
done. On the contrary, strong deviation from the LK be- 
havior is observed for the second harmonics. The main 
feature of this latter component being the zero amplitude 
occurring at a t/b value depending only on the effective 
mass ratio m*^/m*. Finally, it can be remarked that the 
studied ID-chain does not yield frequency combinations, 
only harmonics of the fundamental frequency. In a next 
step, it is planned to consider 2D networks of compen- 
sated orbits which account for the FS of many organic 
metals and are known to give rise to such phenomenon. 
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